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Abstract 

Supersymmetric (SUSY) theories include many new parameters, some of which are CP- 
violating. Assuming that SUSY is found at a future high-energy collider, we examine the 
decay ti — > t\T~r + with an eye to obtaining information about the new CP phases. We 
show that there are two CP-violating asymmetries which can be large in some regions of 
the SUSY parameter space. These involve measuring one or both of the t spins. These 
asymmetries are particularly sensitive to 4>A t > the phase of the trilinear coupling A t . 
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I. INTRODUCTION 



Supersymmetry (SUSY) is usually invoked to solve the hierarchy problem of the 
Standard Model (SM). Because such models are so theoretically compelling, many 
believe that SUSY must be found in nature, and much work has gone into elucidating 
various aspects of such theories. 

In SUSY theories each ordinary fermion and gauge boson has a superpartner, 
respectively of spin and spin |. These models contain a large number of unknown 
parameters. Many of these parameters can be complex, so that, in general, SUSY 
theories predict large CP-violating effects due to the presence of SUSY particles. 
The purpose of this paper is to begin a systematic high-energy exploration of CP 
violation in SUSY couplings. That is, throughout this paper we assume that the 
CP-conserving parameters are known, but that those which violate CP remain to be 
measured. Of course, most measurements will not probe single CP-violating SUSY 
phases. However, if one makes many measurements, then the combined effect will 
be to measure or constrain all SUSY phases. 

The low-energy CP-violating properties of the couplings of superparticles have 
been explored. For example, constraints due to electric dipole moments (EDMs) 
have been studied in detail [1-7]. Contributions to meson mixing have also been 
investigated [8], and the effect of SUSY theories on CP violation in the 5-meson 
system has been examined [9]. 

The EDM experiments are particularly constraining and have led to the so-called 
SUSY CP problem - if sfermion masses are of order the weak scale and complex 
SUSY parameters have phases of order unity, calculations of EDMs generically yield 
values in excess of the experimental limits. Furthermore, the experimental limits 
are expected to be improved substantially over the next few years, which will lead 
to further constraints on the SUSY parameter space [7]. Theoretical calculations 
of EDMs are highly model-dependent, and there are several scenarios in which the 
problem is ameliorated [10]. For example, one simple way to avoid the SUSY CP 
problem is to assume that sfermion masses are of order the weak scale, but that 
all relevant CP-odd phases are highly suppressed. Another is to adopt the opposite 
stance, taking the phases to be of order unity, but sfermion masses to be very large 
(of order several TeV). Scenarios have also been studied in which only the first- 
and second-generation sfermions are very heavy, but the third-generation sfermions 
are near the electroweak scale. In other cases, authors have made assumptions 
concerning universality or alignment of various soft-SUSY-breaking terms. Finally, 
in principle there could be substantial cancellations between different diagrams, 
leading to a reduction in the EDMs. 

Assuming that SUSY particles are found, it will be necessary to measure the 
CP-violating properties of their couplings in order to get a complete picture of 
the physical SUSY model. Several studies have already been performed, showing 
how high-energy observables could be used to measure or constrain various SUSY 
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phases [11, 12]. In the present work we consider CP-odd asymmetries related to a 
particular decay mode of the heavier top squark in SUSY. 

The observables that we consider are sensitive primarily to parameters associ- 
ated with the third-generation squarks. EDMs can in principle constrain quantities 
associated with the third generation. (For example, the two-loop Bar-Zee-type di- 
agrams considered in Ref. [3] can directly constrain the CP-violating parameters 
related to the third generation squarks.) We will assume that the SUSY parameters 
that we use, as well as those not directly involved in our calculation (i.e., those 
corresponding to the first two generations), are such that the various EDM con- 
straints are not violated. Our point of view is that our observables may be used to 
provide an independent measurement of relevant CP-odd SUSY parameters. Such 
measurements would be complementary to the constraints coming from low-energy 
observables such as EDMs or meson mixing. If indeed the EDM constraints end 
up predicting particular values for the parameters we consider, then a high-energy 
measurement should verify such predictions. 

In SUSY theories, there are two Higgs doublets. As a result, there are three 
neutral Higgs bosons - two scalars and one pseudoscalar. If CP is violated, when 
one evolves from the gauge basis to the mass basis, these neutral particles can mix, 
so that the couplings of the physical Higgs bosons are generally mixtures of scalar 
and pseudoscalar. Recently, it was shown that this mixing could be probed in the 
process H° — > ti [12]. The measurement of such a mixing would be one signal of 
CP violation in SUSY theories. The authors of Ref. [12] note that this signal could 
also be seen in H° — > rf if the H° is too light to decay to ti. As such, this method, 
and ours described below, are applicable to high-energy colliders such as the Large 
Hadron Collider or Next Linear Collider. 

SUSY theories also contain two scalar superpartners of the top quark, one for each 
t-quark helicity, called ti and Ir. These two "stops" can mix, resulting in two mass 
eigenstates t 1 and t 2 whose masses can be very different. Here we adopt the standard 
notation > m^. In this paper we examine the decay process t 2 — > Iit~t + . This 
comes about principally through the exchange of an intermediate Z boson or any of 
the Higgs bosons. The measurement of CP violation in this process will therefore 
probe the CP phases in the stop couplings, as well as scalar-pseudoscalar mixing in 
the Higgs sector. 

There are a variety of CP-violating asymmetries in this process, some of which 
depend on strong (CP-conserving) phases. However, since the r's are leptons, they 
cannot emit gluons, and the process t 2 — > tiT~r + does not have QCD-induced strong 
phases (any gluons exchanged between the two stops only serve to renormalize the 
stop couplings). In order to generate strong phases, we therefore include the widths 
of the mediating particles. 

In order to calculate the rate and the various CP-violating asymmetries for t 2 — > 
tiT~r + , various SUSY parameters must be specified: the masses and couplings of m~ tl 
and mi 2l and the masses, mixings and widths of the Higgs bosons. However, these 
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parameters are not all independent - they can be computed from the underlying 
SUSY parameters. To be specific, there are 7 fundamental SUSY parameters which 
strongly affect this decay process (see Sec. IV for an explicit list). 

However, since the aim here is to measure CP violation, we are mainly interested 
in the underlying CP-violating SUSY parameters. We note that all CP-violating 
asymmetries in t 2 — » Iit~t + depend only on the stop couplings, the mixings of the 
Higgses, and their widths. These quantities in turn depend on several CP-violating 
(i.e., complex) parameters in the underlying theory - the various trilinear couplings, 
the gaugino masses and fi. Thus, by specifying the phases of the various complex 
parameters, one can predict all CP-violating effects in t 2 — > t\T~r + (recall that 
we have assumed that all CP-conserving quantities are known). We shall assume 
throughout that ji is real and positive. This may be arranged by rephasing the 
various complex parameters [5]. With this phase convention, the asymmetries in 
t 2 — > tiT~r + have a strong dependence on (the phase of the trilinear coupling 
A t ) and only a much weaker dependence on the other phases. (See Sec. IV for a more 
detailed discussion of this point.) Thus, the measurement of CP-violating effects in 
i 2 — > i\T~T + will allow us to measure/constrain (f>A t - 

A computer program - "CPsuperH" - has been written by others to generate 
the physical SUSY parameters (including the masses, mixings and widths of the 
Higgs bosons, the masses of the stops and the couplings of the Higgs bosons to the 
(s)fermions) from the underlying SUSY parameters [13]. Thus, with this program in 
hand, it is necessary only to specify the values of the fundamental SUSY parameters. 
We use CPsuperH extensively in our numerical work below and have adopted their 
notation for the coupling constants. 

In Sec. II, we compute the amplitude for t 2 — > t\T~T + , as well as its CP-conjugate. 
The various CP-violating signals are examined in Sec. III. We show that there are 
two CP asymmetries which can be large in some regions of SUSY parameter space. 
These involve the measurement of one or both of the r spins. In Sec. IV, we present 
the predictions for these asymmetries. The measurement of these asymmetries will 
provide information about the SUSY CP phases. We conclude in Sec. V. The 
Appendix contains additional information about the various CP asymmetries. 

II. AMPLITUDE FOR t 2 i 1 r-T+ 

We now present the details of the t 2 — > t 1 T~T + calculation. There are box- 
diagram contributions to this process, but they are expected to be small. Thus, 
this decay is mediated principally by the Z boson and by any of the three neutral 
Higgs bosons. That is, we have t 2 — > t\{Z,Hj) followed by (Z,Hi) — > t~t + . In 
the following we always take the t~-t + invariant mass to be above the Z pole, 
and so we assume that only neutral Higgs-exchange diagrams are important. While 
interferences between the Higgs diagrams and the off-shell Z could in principle give 
contributions to the CP asymmetries we consider, these contributions are expected 
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to be small, unless the lightest Higgs boson has a mass very close to that of the Z. 
To simplify our discussion, we ignore any contributions from the Z} 

As mentioned in the introduction, the scalar and pseudoscalar Higgs particles 
in SUSY can mix if CP is violated, leading to three physical neutral Higgs bosons 
that do not have well-defined CP transformation properties. In general, the scalar- 
pseudoscalar mixing need not be small. The squarks also mix, and in general the 
up-type squark mass-squared matrix is 6 x 6. However, it is usually assumed that 
the up-type squarks are block-diagonalized by the same unitary transformations 
that diagonalize the up-type quarks. In this case, one obtains a simplified form for 
the mixing matrix: one only needs to diagonalize a set of 2 x 2 matrices. We adopt 
this simplifying assumption for the stops and consider only the mixing between 
and t R . The CPsuperH program calculates the masses and mixings of both the stops 
and the Higgs bosons. 

The couplings between the stops and Higgs bosons may be described by the 
following Lagrangian: 

£hu = v J2 9Hj*t k Hit*t k , ( 1 ) 

where v is related to the vacuum expectation values of the two Higgs doublets, 
v = \Jv\ + v%, and where % — 1, 2, 3 and j, k = 1,2. The coupling constants gn t t*t k 
are defined as follows [13], 

vg H ^ k = (T^)^O m U^ k , (2) 

where O and IP are the Higgs and stop mixing matrices, respectively. The three 
2x2 matrices r a *** depend on the SUSY parameters A t , fi, cos (3 and sin/5 (where 
tan/5 = V2/V1). Expressions for the matrices T at * 1 may be found in Appendix B 
of Ref. [13]. In the CP-invariant limit, the couplings involving scalar Higgs bosons 
are real and those involving the pseudoscalar Higgs are purely imaginary. If CP is 
broken, these couplings are in general complex. 

We also require the couplings between the tau leptons and the various Higgs 
bosons. The relevant piece of the Lagrangian is 

Chtt = ~9r Yl R i T {dk-T + WHifrl 5 ) T , ( 3 ) 

i 

where % = 1,2,3. The constant g T and the scalar and pseudoscalar coupling con- 
stants, gfj. Tr and g^TTi respectively, are real. At tree-level the couplings are given 



In order that this approximation be reasonable, we will always consider cases in which the lightest 
Higgs is relatively well-separated in mass from the Z. If Z contributions end up being required, 
they could be incorporated in a straightforward manner. 
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by the expressions [13] 



tree 9 m r S.trcc On Ptree r\ , r, i A \ 

9 * =W fe= ^' fe=-0 3 ,tan^. (4) 

The program CPsuperH produces quantum-corrected versions of the above expres- 
sions. In the limit that CP is conserved, a given Higgs boson Hi has either a scalar 
or a pseudoscalar coupling, but not both. 

The final ingredient in our calculation is the Higgs-boson propagator. In the 
decay t 2 — > tiT~T + , we allow one or more of the Higgs bosons to go on-shell. The 
reason for this is two-fold. First, this will dramatically increase the rate, giving a 
possibility that experimentalists will in fact be able to measure the asymmetries that 
we construct. Second, as noted above, the strong phases in this process are obtained 
mathematically by including the widths, i.e. absorptive pieces in the propagator. 
These absorptive parts have a much larger effect near resonance. 

There is one complication here. It is true that CPsuperH computes all the physical 
properties of the Higgs bosons. In particular, the width of the Hi boson is found by 
calculating the 1-loop contributions of the light particles (taken here to be principally 
b and r), and applying Cutkowsky rules. However, in our calculation the "off- 
diagonal widths" come into play. That is, this 1-loop diagram can also be used 
to relate a Hi boson to a Hj boson, and the absorptive parts of such diagrams 
can be important for those CP asymmetries in our decay process which rely on 
strong phases. CPsuperH does not give these off-diagonal widths, but they have been 
computed in Ref. [14]. 2 A good approximation for the Higgs propagator (divided 
by T) is 

D(M 2 ) = 

/ M 2 — m 2 Hl + i Imfin ilmn 12 ilmR 13 \ _1 

% Irnfisi M 2 - m\ 2 + % Imfl 22 % Imn 23 , (5) 

V i Imn 31 i Imn 32 M 2 - m 2 Hz + i lmU 33 J 

where M 2 = {p 1 + p 2 ) 2 , with p 1 and p 2 being the four-momenta of the r" and t + , 
respectively. Explicit expressions for the absorptive parts of the Higgs boson self- 
energies, Iinn^M 2 ), may be found in Ref. [14]. For our purposes, it is sufficient to 
consider only the contributions to ImIljj(M 2 ) arising from the b-b and loops 3 , so 



2 CPsuperH does include off-diagonal contributions to the Higgs self-energy in its calculation of 
the Higgs boson masses [15]. 

3 We restrict our consideration to cases in which the neutral Higgs bosons can only decay into 
(non-i) quark and lepton pairs. In principle, M 2 could still be large enough that other terms 
should be included in the calculation of Imlljj (M 2 ). In practice, however, we are only interested 
in the effect of ImUij(M 2 ) near resonances, i.e., when M 2 ~ m 2 Hi . In such cases it is sufficient 
to restrict our attention to the b and r contributions. 
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FIG. 1: Feynman diagram for the decay t 2 —> t\T~T + . Note that the Higgs propagator has 
off-diagonal terms, so transitions Hi — > Hj are allowed. The shaded circle at the Hj-r-T 
vertex denotes the fact that the effective coupling calculated by CPsuperH includes some 
loop corrections. 



that 

where 

Imng(M 2 



imn- 



(M 2 ) 




+ imn- 



(M 2 ), 



(6) 



p p 
9hm>9hM 



AV 2 (l,« 6 ,« 6 ), (7) 



4«r) 9%- T 9%- T + 9Herr9 P H,fr ^ «r, «r) , (§) 



with A(x, y,z) = x 

S,P 



y 2 + z 2 — 2(xy + xz + yz) and kj, jT = m 2 T /M 2 . The constants 
(ft, and S 1 ^, are associated with the Hjbb couplings and are denned in analogy 
with Eq. (3). (The tree-level expressions for the coupling constants are identical 
to those given in Eq. (4), but with the substitution m T — > m b .) Note that the 
diagonal elements are related to the widths of the Higgs bosons through the optical 
theorem, ImIIjj(m^.) ~ r(iJj)m^, giving rise to the usual Breit-Wigner form of the 
propagator if one takes M 2 = m 2 H . in the self-energies and ignores the off-diagonal 
terms. We include the off-diagonal terms in our calculation, since in some cases they 
have a non-negligible effect on the asymmetries that we consider. 

Considering only the contributions due to the Higgs diagrams (see Fig. 1), the 
amplitude for the decay £J — > t^T~T + is given by 

A(si, s 2 ) = Bu(pi, si)v(p 2 , s 2 ) + Cu(p 1 , s 1 )-f 5 v(p 2 , s 2 ), (9) 

where the indices '1' and '2' refer to r~ and r + , respectively. The complex coeffi- 
cients B and C are given by 

B = vg^QH&hDiAM^gS , (10) 
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C = ivgrYsdH^hDijiM^gP 



Calculation of the absolute value squared of the amplitude, for specific spin states 
si and s 2 , yields the expression, 



l-4(s„s 2 )f 



- \B 

2 1 

+-|Cf 



(M 2 - 4m 2 ) (1 - si • s 2 ) + 2 pi • s 2 P 2 • si 

M 2 (1 + Si • s 2 ) - 2 pi • s 2 • si 
-m T (5C* + 5*C) (pj • s 2 + p 2 • si) 



(12) 



in which we have adopted the convention eoi 2 3 = +1. 

In order to calculate CP- violating asymmetries we must determine the amplitude 
for the CP-conjugate process. This is 



where 



A = Bu(p!, si)f (p 2 , s 2 ) + Cw(pi, si)7 5 w(p 2 , s 2 ), 



(13) 

(14) 
(15) 



(recall that g T and ^§^ T are real). As one might expect, the CP-conjugate amplitude 
is obtained from the original amplitude by complex-conjugating the weak phases and 
leaving the strong phases unchanged (the strong phases appear in the propagator 



matrix) 4 . The expression for 



A 



is obtained from Eq. (12) by the replacement 

(B,C) — > (B,C). (Note that the explicit 'i' appearing in Eq. (12) does not change 
sign under CP conjugation.) 



4 Actually the situation is somewhat subtle. We have derived the expression for A by determining 
the CP-conjugates of the initial and final states and then calculating the corresponding ampli- 
tude. Alternatively one could reverse the signs on all weak phases <f>i (<pA t , 0A b , etc.) and 
recalculate A. Under the transformation <pi — > —(pi, one finds that zimlli3 — > — iimlli3 and 
ilmn 2 3 — > — ilmn 2 3; i.e., some of the strong phases change sign. As a result, the 1-3 and 2-3 
elements of the Higgs propagator also change sign. A careful consideration of sign changes that 
occur simultaneously in some of the coupling constants {gu 3 Qti ~~ * ~9H 3 t*t 1 ' 9h 3 tt ~~ * ~9h 3 tt> 
g^ iTr — > — g^TT an d 9h 2 ~t ~~ * ~9h 2 tt) snows that the expressions given for B and C are 
precisely correct. 
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III. CP ASYMMETRIES 



We can form three CP- violating asymmetries from the expressions for \A\ and 

2 

A , using various combinations of polarizations to extract different terms. The 
irst asymmetry is the usual rate asymmetry, for which it is assumed that spins are 
not measured. The second asymmetry is a "single-spin" asymmetry. It requires the 
measurement of a single r spin in the decay and is sensitive to the real part of BC*. 
The third asymmetry is a triple-product asymmetry and requires the measurement 
of both spins. It is sensitive to the imaginary part of BC* . 

In general, the construction of the CP asymmetries proceeds in two steps. In 
the first step, we sum over the r ± spins (possibly in an asymmetric manner) for the 
process — > ^[t~t + . We denote the resulting expression for the amplitude squared 
by U \A\ 2 ," where '(a)' refers to the specific prescription employed when summing 

(a) 

2 



over the spins. In the second step we construct 



A 



the analogous quantity 

(a) 

for the CP-conjugate of the process, and subtract it from \A\ 2 ^ . The resulting 

expression is odd under CP. We also integrate over the invariant mass M of the r ± 
pair, starting at some point above the Z resonance but including one or more Higgs 
resonances. The CP asymmetries can in principle be increased by a judicious choice 
for the range of this integration. Assuming M > M mm , where M min > 2m T , 5 we 
generically define the width associated with the a th prescription for summing over 
spins as 



(a) = 1 [ m h ~ m h dM f 2 

^min " 1287T3m? J V in M V ' 



(a) 



X 



t 2 " J "min 

1/2 



\(m{,ml,M 2 )\(M 2 ,m 2 T ,m 2 T )} . (16) 



The width for the CP-conjugate process is obtained by the replacement (B, C) — > 
(B, C) and is denoted Y^f! . 

\ i 1M min 

A. Rate Asymmetry 

The first CP asymmetry is the usual rate asymmetry and is obtained by summing 
symmetrically over the spin states s\ and s 2 to yield 

|^| 2 = ]T \A(s 1 ,s 2 )\ 2 = 2\B\ 2 (M 2 -Am 2 T )+2\C\ 2 M 2 . (17) 

rate V / 

spins 



5 In practice, we always assume that M mm is above the Z resonance so that we can ignore Z 
contributions. 
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This expression may be inserted in Eq. (16) to obtain r^ min , and similarly for ^M m[n 
(we drop the 'rate' superscripts in this case). The resulting CP asymmetry is then 

r — r 

/irate _ ^min -^min (18) 

CP T1 , Tl V / 

M min + M min 

A non-zero value for this asymmetry requires both a weak-phase difference and 
a strong-phase difference between the diagrams contributing to the decay. Since 
we obtain our strong phases from absorptive parts of the Higgs boson propagator, 
simultaneously non-negligible values for the asymmetry and the rate itself require 
the interference of two or three nearby Higgs boson resonances. 

We have found numerically that the rate asymmetry tends to be very small for 
t 2 — > tiT~T + , with values typically at the sub-one percent level. There are a few 
factors that conspire to make the rate asymmetries small; these are described in the 
Appendix. At this point we simply note that this particular observable will not be 
a very useful tool for the study of SUSY CP violation in this decay process. It turns 
out, however, that the measurement of the spin of one or both final-state leptons can 
lead to a dramatic increase in the sensitivity of t 2 — > tiT~T + to CP phases. The two 
asymmetries described below use lepton spins to probe the cross-terms containing 
B and C in Eq. (12). 6 These asymmetries can be quite large numerically. 



B. Single-spin Asymmetry 



We first construct a CP asymmetry that uses the polarization of only one of the 
leptons. This single-spin asymmetry extracts the piece in that is proportional 
to the real part of BC* . In constructing this asymmetry, it is useful to calculate 
|-4( s i) s 2 ) | 2 explicitly for the different spin combinations of the t ± . The spin four- 
vectors Si and s 2 are determined in the usual manner by Lorentz-transforming unit 
three-vectors from the rest frames of the particle and antiparticle. Let the unit 
three- vector s denote a specific spin state in the rest frame of the r + or r~ . Lorentz- 
transforming s to a reference frame in which the r + or r~ has four-momentum 
p M = (E,p) yields 

V m m(E + m) ) 



We present our results as asymmetries depending on the spins of the r~ and r+. In practice, 
experimentalists would use the decay products of the taus to probe their spins. In that case, 
one could construct CP asymmetries depending on the momenta of the decay products. See 
Rcf. [16] for an explicit comparison of 'theoretical' triple-product asymmetries involving the 
D* polarization states in B — > D*£U with more 'experimental' asymmetries that use only the 
momenta in the four-body final state of the decay B — > (D* — > Dir)tD. 
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For a given set of momenta in the process t 2 — > tir~r + there are four spin combi- 
nations that one could consider, corresponding to the four combinations of ±si and 
±s 2 - In constructing the single-spin asymmetry it is convenient to work in the helic- 
ity basis, taking the unit three-vectors to be aligned or anti-aligned with the lepton 
three-momenta in the rest frame of the lepton pair. This choice also maximizes the 
single-spin asymmetry. (In the following subsection we will construct triple products 
using the spins si and s 2 . In that case it will be convenient to employ a basis in 
which the unit three- vectors are perpendicular to the lepton three-momenta.) 

Since \ A\ 2 is Lorentz- invariant, we may choose any convenient frame of reference. 
We work in the rest frame of the r ± pair, so that p*_^ 2 — {B, with p denoting 
the three-momentum of the r~ in that frame 7 . Defining s\ = p/\p\ and Lorentz- 
transforming ±§i from the r~ rest frame to the Higgs rest frame, we have Si ±M = 
± (|p| /m T , Ep/(m T \p\)). Similarly, for r + spins that are in the ±p/\p\ direction in 
the r + rest frame, we have s 2 ±M — ± (— \p\ /m T , Ep/(m T \p\)) in the Higgs rest frame. 
Explicit calculations for the four spin combinations yield 

|-4(±,±)| 2 = 0, (20) 
\A(±,T)\ 2 = \B\ 2 (M 2 -4m 2 T ) + \C\ 2 M 2 



=F (BC* + B*C) My M 2 - Am 2 T ■ ( 21 ) 

Summing all four spin combinations reproduces the result in Eq. (17), as expected. 
Taking the difference of the two non-zero spin combinations yields 



l^l 2 



spin 



|^(+, -)| 2 - |^(-, +)| 2 = -4 Re (BC*) M^/M 2 - 4m* . (22) 



The resulting CP asymmetry is given by 

■pspin j=f s P m 

^spm _ M min M min ^23) 

^ M ■ r 'm ■ 

1 ™min J "min 

in which we have normalized the asymmetry using the total widths V M . and V M . . 
Note that this asymmetry only requires the measurement of a single spin, since the 

++ and combinations are zero. As was the case for the rate asymmetry, this 

single-spin asymmetry requires the interference of two amplitudes with a non-zero 
relative strong phase. 

In contrast to the rate asymmetry, the single-spin asymmetry can be fairly large 
numerically since it does not suffer from the same cancellations as those of the rate 
asymmetry. 



7 In the remainder of this work we will refer to the rest frame of the pair as the rest frame of 
the decaying (virtual) Higgs. 
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FIG. 2: Definitions of angular variables and unit vectors used in the calculation of the 
triple-product CP asymmetry. The unit vector z specifies the direction of the virtual Higgs 
as seen in the ti rest frame. The vector p = pn specifies the r~ momentum in the Higgs 
rest frame. 



C. Triple-product Asymmetry 

For the triple-product CP asymmetry we must use a different basis for the spins 
si,2- We work in the rest frame of the virtual Higgs (i.e., in the rest frame of the 

pair), denning the +z direction in that frame to be the direction of the Higgs as 
seen in the ti rest frame. We let the unit vector z denote this direction (see Fig. 2). 
We let the unit vector n denote the direction of the r~ in the Higgs rest frame. 
To specify the spin four- vectors of the t^, we start by denning the corresponding 
unit three- vectors in the individual rest frames of the t~ and the r + . These are 
then Lorentz-transformed to the Higgs rest frame. In the individual t~ and r + rest 
frames we define the '+' and '— ' spin states as 

. z — n cos 9 ,z — ncos9 

sf = ±— = ± — — , (24) 

\z — ncos9\ smv 



= ±7^r = ±^ . (25) 
\z x n\ sin 9 



where 9 is defined in Fig. 2. Then the vectors sf, sf and n are mutually per- 
pendicular. The 'e' term may then be extracted from Eq. (12) by the following 
prescription: 



i^i 2 | TP = -^J M+> +)i 2 + -)i 2 - -)i 2 - +)r 



= Aim {BC*) M\j M 2 - Am 2 T , (26) 
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where the angular integration is performed in the Higgs rest frame, and where the 
± spins are defined as in Eqs. (24) and (25). (This basis is different than that 
employed in Eqs. (20) and (21).) Inserting Eq. (26) in Eq. (16), and similarly for 
the CP-conjugate process, we form the CP-odd triple product asymmetry as follows: 

■pTP _ p TP 

^TP = M min _M min , . 

CP -p i -p v / 

1 M min + 1 M min 

2 

Note that the prescription used to define \A\ Tp is equivalent, in this spin basis, to 
the definition 

r AVm = rM min ( 6 (Pl' S l'^, S 2 ) > 0) - T M miQ (e(pi, S U p 2 , S 2 ) < 0) . (28) 

The triple-product asymmetry does not require a relative strong phase between 
interfering amplitudes. In fact, the triple-product asymmetry does not even nec- 
essarily require the interference of two different Higgs resonances - one resonance 
suffices as long as there is scalar-pseudoscalar mixing among the Higgs bosons. 



IV. RESULTS 

As has been noted above and discussed in the Appendix, the rate asymmetry in 
t 2 — > t\T~r + will in general be very small, so we do not consider it further here. 
The single-spin CP asymmetry can be large, but requires the interference of two or 
three nearby Higgs resonances. If the resonances are spread far apart, the product of 
any two interfering amplitudes will be suppressed, since one amplitude will be much 
smaller than the other. The strong-phase difference in this CP asymmetry is due to 
the width difference between the contributing Higgs bosons as well as the relative 
values of their masses (see Eq. (A. 7)). The triple-product CP asymmetry does not 
require a strong-phase difference and can receive large contributions from single 
resonances (due to scalar-pseudoscalar mixing) and from the interference between 
two resonances. The Appendix contains further discussions of these two distinct 
types of contributions. 

In SUSY it is somewhat natural to have nearby Higgs resonances, so that it is 
not unreasonable to hope for large single-spin and triple-product asymmetries. At 
tree level there is a well-known relation between the masses of the heavier scalar 
Higgs (H) and the pseudoscalar Higgs (A), 

m 2 H -m 2 A ~ m| sin 2 2(3 , (29) 

a relation which is valid in the limit m A 3> m| [17]. For large tan/5, sin 2 2/3 
and thus m 2 H m 2 A . If CP is not conserved, loop corrections lead to the mixing of 
the scalar and pseudoscalar Higgs bosons [18]. This mixing can lead to the breaking 
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of the near degeneracy in masses, as has been explored in some detail in Ref. [19]. 
Nevertheless, for some combinations of parameters (in particular, it appears, if \fx\ 
is not too large) there remains a near degeneracy of the two heavier Higgs bosons. 

In the following we choose a few points in the SUSY parameter space and use these 
to illustrate some of the features of the single-spin and triple-product asymmetries. 
Our treatment is certainly not exhaustive, but is meant simply to point to the 
usefulness of using these asymmetries as a tool for exploring SUSY CP violation. 
We choose values for the charged Higgs mass and tan (3 that are consistent with the 
recent bound from Belle, tan/5 /m H ± < 0.146 GeV -1 [20]. 

The asymmetries A S q P u an d ^cp defined in Eqs. (23) and (27) have been inte- 



grated over M, where M — y(pi + P2) 2 is the invariant mass of the pair (see 
Eq. (16)). Since the r's are produced by the s-channel decays of Higgs bosons (see 
Fig. 1) , the differential widths (in M) for the decays tf — > tfr~r + are expected 
to exhibit resonant peaks in the vicinities of the neutral Higgs masses; i.e., when 
M m rriH v with % = 1,2,3. The CP asymmetries that we have constructed have 
a sensitive dependence on the resonance structure of the decay. For example, as 
mentioned above, the single-spin asymmetry requires the interference of different 
resonances, and will be suppressed if the resonances are too far apart. The triple- 
product asymmetry can receive contributions from single resonances, but in some 
cases the contributions due to different resonances could be of opposite signs, lead- 
ing to a suppression of the (integrated) quantity A^p. Experimentalists may wish 
to enforce various cuts on the integration over M in order to optimize both the CP 
asymmetry and its statistical significance. 

To display some of the above features more clearly, let us define several differential 
quantities. First, we define dT SUTa /dM to be the sum of the total (i.e., summed over 
spins) differential widths for the process and the anti-process. That is, 



A plot of dT snm /dM as a function of M will show resonance structure, with peaks 
appearing at the Higgs masses. Furthermore, integrating dT sum /dM over M gives 
the sum of the widths for the process and the anti-process, 



r sU m is the quantity that appears in the denominator of the CP asymmetries defined 
in Eqs. (23) and (27). Recall as well that it has an implicit dependence on the lower 
limit, M min , of the integation in the above expression. We also define two CP- 





(30) 




(31) 
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FIG. 3: Plot of ( solid line ) 5 ( long dashed line) and (short dashed line) 

for a case in which the two heaviest Higgs resonances (near 126.1 GeV and 128.5 GeV) 
have a significant overlap. 



asymmetric differential widths, 



dY 



(a) 
CP 



dM 



128vr 3 m? M 
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(a) 



A 



(a), 



X 



A (m|,m|,M 2 ) \(M\m 2 T ,m; 



1/2 



(32) 



where 'a' corresponds to 'spin' or 'TP'. When integrated over M, these differential 
widths form the numerators in the expressions for the CP asymmetries in Eqs. (23) 
and (27), so that 



A 



(a) 
CP 



-f 

sum JM min 



^ 2 - m h dT ( cl 
dM 



dM . 



(33) 



In the following we will show plots of dT^/dM, dT^/dM and dT snm /dM. The 
relation between these plots and the integrated CP asymmetries and A^p 

follows from Eq. (33): The spin (TP) asymmetry is the ratio of the area under the 
curve dr S Q P n< " TF ^/dM to the area under the curve dT sum /dM. 
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Figure 3 shows plots of the differential widths dT S Q F n /dM, dTj^/dM and 
dT snm /dM for a particular set of input parameters. To generate these plots we 
have used CPsuperH to determine the masses and couplings of the Higgs bosons 
as well as the masses of the stops. Some of the input parameters are as follows: 
m H ± = 150 GeV, fj, = +200 GeV, \A t \ = 600 GeV, <f) At = 90°, tan/3 = 20, 
m,Q 3 = 350 GeV and m^ 3 = 400 GeV, where we have followed the notation of 
CPsuperH. (The decay process strongly depends on all of these parameters. Other 
parameters have been chosen such that the neutral Higgs bosons cannot decay into 
SUSY particles.) The only non-zero phase is that of the trilinear coupling A t . The 
masses of the stops are determined by CPsuperH to be ~ 259 GeV and ~ 
511 GeV. Similarly, the neutral Higgs masses (widths) are mu x — 114.38 (0.15) GeV, 
m H2 ~ 126.12 (1.09) GeV and m H . A ~ 128.45 (0.95) GeV. 

As noted above, to obtain the integrated CP asymmetries, we determine the areas 
under the curves in Fig. 3 and calculate the appropriate ratios. A comparison of the 
area under the dashed curve (dT^p/dM) to that under the solid curve (dT sum /dM) 
indicates that the CP asymmetry will be relatively large in this case. A similar 
analysis for the single-spin asymmetry leads one to the conclusion that that asym- 
metry will be somewhat smaller. Our numerical results confirm these qualitative 
observations. Integrating over all three Higgs resonances we have A^p 1 ~ —14% 
and Al? ~ -31%. 

Consider also a few other technical details associated with this example: 

1. One might wonder why the resonance evident in dT sum /dM near the heaviest 
Higgs at 128 GeV is larger than that near 126 GeV. It turns out that the 
difference in the strengths of these resonances is due primarily to differences 
in the Higgs-stop-stop couplings. 

2. In this example the two heaviest Higgs resonances have a relatively significant 
overlap, allowing for some amount of interference between the resonances. 
The single-spin asymmetry requires such interference, while the triple-product 
asymmetry can also receive contributions from single resonances. 

3. Including only the heavier two resonances in the integration over M increases 
A S Q P n (in terms of its magnitude) to about —19%, but decreases A^p (again 
in terms of its magnitude) to about —24%. The increase in the single-spin 
asymmetry is not difficult to understand: The numerator in Eq. (33) is nearly 
unchanged (because the Hi resonance is so far separated from the other two 
and does interfere very much), but the denominator decreases from r sum ~ 
0.13 GeV to r sum ~ 0.098 GeV. 

4. Observation of a triple-product asymmetry in the vicinity of a single, widely- 
separated resonance (such as Hi in this case) would be evidence of scalar- 
pseudoscalar mixing in the Higgs sector. 
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FIG. 4: Plot of ^ m (solid line), (long dashed line) and £fi (short dashed line) 

for a case in which the two heaviest Higgs resonances are somewhat spread apart. The 
input parameters for this plot are noted in the text. 



A second example is shown in Fig. 4. The input parameters for this example 
are the same as those of the previous example, except that ji and A t have been 
increased in magnitude (/i = +500 GeV and \A t \ = 700 GeV in this case). As 
a result, and m# 3 have moved farther away from each other and mH 1 has 

increased. (The neutral Higgs masses (widths), as determined by CPsuperH, are 
m Hl ^ H6.09 (0.16) GeV, m H2 ~ 124.73 (0.83) GeV and m Hz ^ 130.11 (0.74) GeV, 
in this case. The stop masses are ~ 228 GeV and m t ~ 2 ~ 526 GeV.) The CP 
asymmetries and integrated width in this example, when integrated over all three 
(only the heavier two) resonances, are Aq^ 1 ~ —5% (—6%), Aqp ~ —36% (—30%) 
and r sum ~ 0.93 (0.83) GeV. The single-spin asymmetry has been reduced some- 
what, compared to the previous example, due in part to the fact that the resonances 
are more widely separated in this case. Also note that the triple-product asymmetry 
suffers from some amount of cancellation between the second and third resonances, 
since the area under the curve near these two resonances is of opposite sign. In such 
a scenario it might be experimentally preferable when measuring the triple-product 
asymmetry to restrict the integration over M to smaller regions or to weight the 
region near the middle resonance with a relative negative sign (to avoid the cancel- 
lation) . 
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FIG. 5: Plot of the integrated triple-product and single-spin asymmetries, A^p and A S Q P n , 
respectively, as functions of 4>A t - The input parameters are the same as those used for the 
plot in Fig. 3, except that <pA t is allowed to vary. The solid (dashed) curve corresponds to 
t4^p integrated over all three (only the heavier two) resonances. Similarly, the long-short- 
dashed (dashed-double-dotted) curve corresponds to A s ^ F n integrated over all three (only 
the heavier two) resonances. 



Figure 5 shows the variation of the integrated CP asymmetries as functions of 
the phase of A t , <j>A t , for the same set of input parameters used in Fig. 3 (except, of 
course, for 4>A t )- Note that the masses of the Higgs bosons and stops, as well as the 
widths of the Higgs bosons, are themselves functions of (f>A t , so their values change 
over the course of the plot. The four curves correspond to A^p and A S Q P n , with two 
different choices for M m j n , one of which includes all three resonances and one of which 
includes only those corresponding to H 2 and H 3 . In this example, including the first 
resonance tends to decrease the magnitude of A S Q P n , while it sometimes increases the 
magnitude of A^p and sometimes decreases it. This is evident from the plot. Values 
as large as about 64% and 49% are observed for the triple-product and single- 
spin asymmetries, respectively. It is possible to obtain even larger triple-product 
asymmetries than those shown in Fig. 5. For example, using the same parameters 
as those used for Fig. 4, but taking <pA t — 135°, one obtains A^p ~ —92% and 
r sum ~ 0.73 GeV when integrating over all three resonances. 

We have performed a limited study of the effects of allowing other SUSY pa- 
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rameters to have non-zero phases. The phase in the trilinear coupling Aj, could in 
principle affect our asymmetries through the mixing of the Higgs bosons. In prac- 
tice, we have found that this phase only seems to have a non-negligible effect when 
\Ab\ is quite large - of order 5-10 TeV in the cases studied - due to the difference in 
the top and bottom Yukawa couplings. The phases of the gaugino masses could also 
in principle give contributions due to their effects on Higgs-fermion-fermion effective 
vertices. Taking the magnitudes of M 1; M 2 and M 3 to be approximately 200, 400 
and 1000 GeV, we have found that in the cases considered these quantities can make 
changes to the triple product asymmetry that are of order ±0.1. The effect on the 
single-spin asymmetry tends to be smaller in magnitude. 

We have also considered scenarios involving different values for tan/3 or A t . The 
triple-product and single-spin asymmetries are strongly dependent on the particular 
values chosen for these quantities. In the few examples we studied, decreasing tan f3 
from 20 (the value used in Figs. 3-5) tended to decrease the asymmetries. We 
also increased tan/9 in one case (and increased m#± to respect the experimental 
bound on tan j3/m H ±) and found that the triple- product asymmetry could increase 
in that case. Varying the magnitude of A t (while holding its phase constant) can 
also affect the asymmetries. In the cases studied, the asymmetries tended to be 
somewhat constant as \A t \ was varied between 400 GeV and 700 GeV. There were 
more significant changes in the range 700 GeV to 1 TeV. 

We should also make a comment regarding the invariance of our results under 
rephasing of the complex SUSY parameters. As is well-known, it is possible to rotate 
away two phases associated with the various complex parameters in supersymmetry 
(see, for example, Ref. [5]). CPsuperH appears to use one of these degrees of freedom 
to force the SUSY parameter m\ 2 to be real [18, 21]. The remaining complex 
quantities within CPsuperH are the gaugino masses (Mj, with i = 1,2,3), the third- 
generation A-terms (A a , with a = b,t,r) and /i. Six physical CP phases may 
be constructed from these seven complex parameters. These may be taken to be 
Arg(AlMi) and Arg(fiA a ) [5]. (Although our scenario is not identical to that in 
Ref. [5], we may borrow the results from that reference because for our calculation 
we may consider the first two generations to be decoupled from the third. The more 
general case is considered in Ref. [22].) In this work we have used the remaining 
available phase rotation to set Arg (/z) = 0, in which case the phases of the A-terms 
are physical. Nevertheless, an important check of our asymmetries would be to 
study their behavior under a rephasing transformation. For example, if one wished 
to make A t real and positive through the transformation A t — > A t e~ l ^ A t, then the 
phase <pA t would appear in the other complex parameters, since one would also need 
to take fi — > /ie^ A t , A bjT — >■ A biT e~ i<l>A t and Mj — > Mje~^ A * . The asymmetries should 
be unchanged under such transformations. We have performed a numerical check of 
the rephasing invariance of our asymmetries in a few particular examples and have 
found them to be invariant to within the expected degree of numerical precision. 

To summarize the numerical work, we note that while the regular rate asymmetry 
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is expected to be very small for t 2 —> t l t~t + , wc have seen that the single-spin and 
triple-product asymmetries can be quite large. In this work we have seen examples 
where the former can be of order 50% and the latter of order 90%. Such large values 
indicate that these asymmetries provide a promising avenue for the measurement of 
CP violation within SUSY. In particular, the observables studied here are very sen- 
sitive to the phase of the trilinear coupling, A t , as is evidenced in Fig. 5. We have 
also seen that the differential CP asymmetries can have an interesting functional 
dependence on M, the invariant mass of the r 1 * 1 pair. Depending on the particu- 
lar set of circumstances, experimentalists may wish to consider certain resonances 
separately, or to find other ways to optimize the CP-odd signals without sacrificing 
their statistical significance. 

V. CONCLUSIONS 

Many physicists believe that supersymmetry (SUSY) will be discovered at future 
high-energy colliders. If so, first measurements will involve CP-conserving quantities 
such as the masses of SUSY particles, etc. However, SUSY also contains a number 
of CP-violating parameters. It is only through their measurement that one will be 
able to identify the type of SUSY theory which is found in nature. 

In this paper we have studied CP violation in the decay t 2 — > 1\t~t+ . We have 
shown that two CP asymmetries can be quite large in some regions of the SUSY 
parameter space (up to 90% for the parameters we have chosen). They are the 
single-spin and triple-product asymmetries, and involve the measurement of one or 
both of the r spins. Both of these asymmetries depend almost entirely on a single 
CP- violating parameter, the phase of A t , <pA t - Thus, the measurement of the CP 
asymmetries in this decay will allow one to extract or constrain 4>A t - (Future work 
will involve other channels, sensitive to other SUSY phases.) 
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APPENDIX: SOME APPROXIMATE EXPRESSIONS 

In this appendix we provide some approximate expressions that help to clarify 
the discussions of the three asymmetries. 
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1. Rate Asymmetry 



In the text it is noted that the rate asymmetry is suppressed for t 2 — > tiT~T + . 
To see why this is the case, consider the difference in the amplitudes-squared (which 
appears in the numerator of the expression for the rate asymmetry, Eq. (18)), 
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(A.l) 



To derive this expression we have made use of the definition of Imfi^ (M 2 ) in Eq. (8); 
we have also set 



on = arg(^ it -* t - 1 ) . 



(A.2) 



The above expression may be simplified by assuming that the off-diagonal terms 
in the propagator are small (although not negligible) compared to the diagonal 
elements. Expanding Eq. (A.l) under this assumption, we find 
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X lfc = (M 2 - m^. + i Imfi^) (M 2 - m 2 ^ + i lmU kk ) . 
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(A.4) 
(A.5) 



In writing down these expressions, we have omitted a class of terms that could be 
important if all three resonances are close together. 

The expression in Eq. (A.4) is small for a few reasons. First, a cancellation has 
occurred so that only self-energy loops with b quarks contribute to the required 
strong phase difference (the r contributions disappear in going from Eq. (A.3) to 
Eq. (A.4)) 8 . Second, the terms Imn^ImlT^ and — Imll^Imll^ tend to cancel each 

S P S P 

other. Part of the reason for this is that the coupling constants ^ and g#.- r are 
equal at tree-level. 



Note the importance of keeping the off-diagonal propagator terms in this case. 
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2. Single-spin Asymmetry 



In our numerical work we keep the full 3x3 Higgs propagator, but to understand 
some of the dynamics for the single-spin asymmetry it is sufficient to keep only the 
diagonal terms. In that case we obtain 
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where we have also made the approximation ImILjj(M 2 ) ~ YjUiu^ with Fj = T(Hj 
If we define the strong phases Sj(M 2 ) as follows 
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we have 
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In this form we can see clearly that the single-spin asymmetry depends on the sine 
of the relative strong phase between the interfering resonances. But it is perhaps 
somewhat counterintuitive that the single-spin asymmetry depends on the cosine of 
«fc — otj, the weak phase difference. 

To understand the presence of 'cos(o; fc — <x,)' in Eq. (A. 8), consider first the CP- 
invariant limit. In this case there is no scalar-pseudoscalar mixing - two of the Higgs 
bosons are pure scalars and one is a pure pseudoscalar. Let the indices js correspond 
to the two scalar Higgs bosons and let jp correspond to the pseudoscalar Higgs. Then 
a js = or 7r and otj p = ±7r/2 (i.e., the stops have real (imaginary) couplings to the 
scalars (pseudoscalar) - see Eq. (A.2)). Furthermore, 9 S u ip - T = 9h js tt = (i- e -> the 
Higgs bosons have pure scalar or pseudoscalar couplings to the taus). It is clear, 
then, that Eq. (A. 8) is zero in the CP-invariant limit. Either the product of the 
scalar and pseudoscalar couplings is zero, or the weak-phase difference is ±7r/2, so 
that the cosine is zero. 
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Now suppose we allow CP to be broken by a small amount, so that there are 
two "mostly scalar" and one "mostly pseudoscalar" Higgs bosons. In this case there 
are two distinct types of contributions to the asymmetry. One contribution comes 
from the interference of a scalar and a pseudoscalar, denoted by the indices k and j, 
respectively. In this case we can redefine a,j = for the (mostly) pseudoscalar 

Higgs, so that a.j measures the departure of atj from the CP-invariant limit of ±7r/2. 
Then cos(o;fc — aij) = ±sin(o; fc — dij), and the expected sine of the weak phase 
difference appears. Alternatively, suppose j and k are both mostly scalars, but 
gH.- T is not quite zero. In this case there can be a contribution to the asymmetry 
even if = otj. 

To summarize, the single-spin CP asymmetry can receive contributions both 
from a weak phase difference in the interfering amplitudes, as well as from scalar- 
pseudoscalar mixing of the Higgs bosons. In general there will be some combination 
of these effects. Both types of contributions require a strong phase difference. 



3. Triple-product Asymmetry 



Making the same approximations for the Higgs propagator as was made in the 
previous subsection, we find the following expression in the case of the triple-product 
asymmetry, 
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(A.9) 



This expression bears some similarity to that derived for the single-spin asymmetry 
(Eq. (A.8)), but with two important differences. In the first place, the cosine of 
the relative strong phase appears - this asymmetry does not require a strong phase 
difference. In the second place, there is a new term compared to the single-spin 
asymmetry - the term with a single sum over k. Since there is no need to have 
a relative strong phase for the triple-product asymmetry, it is possible to receive 
contributions from single resonances. Such a contribution contains an implicit de- 
pendence on CP-violating phases because it is only non-zero if gfi k - T and Tt are 
both non-zero; i.e., if the Higgs boson in question has both scalar and pseudoscalar 
couplings. This can only happen if CP has been broken. Note that this first term is 
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somewhat analogous to the asymmetry considered by Valencia and Wang in Ref. [12]. 
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